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THE  GAMMA  AND  BETA  FUNCTIONS 
By  W . Edwards  Demi ng 

Definition  of  the  Gamma  function."-  The  Gamma  func- 
tion for  positive  values  of  n can  be  defined  by  the  inte- 
gral 1 

r(n)  = /x"--*  e~x  dx.  (l) 


T(n)  is  read  "the  Gamma  function  of  n. " It  is  truly  a 
function  of  n,  that  is,  it  depends  on  n for  its  value. 

With  the  upper  limit  infinite,  as  here  written,  the  Gamma 
functions  is  said  to  be  complete , but  this  adjective  is 
usually  omitted  except  when  it  is  desired* to  distinguish 
a complete  from  an  incomplete  Gamma  function.  The  incom- 
plete  Gamma  function  is  defined  "by2 

<*>  r (n)  = /Xx"-Ie-r  Ax,  (2) 

CS  X o 

uu 

■ • 

which  is  the  same  integral  except  that  the  upper  limit  is 
a variable.  The  incomplete  Gamma  function  Tx (n)  is  thus 


The  integral  given  here  ae  the  definition  of  the  Gamma  function 

was  studied  by  Euler,  mainly  in  the  form  f*( In  l/x)  ^ dx  (vol.  k of 

o 

his  Inatitutlonis  Calculi  Integralisj  Petrograd,  1770) , and  is  some- 
times referred  to  as  "Euler's  second  integral."  An  excellent  presen- 
tation of  the  Euler  integrals  is  given  by  Joseph  Edvards  in  his  Inte- 
gral Calculus,  vol.  2,  Ch.  2k  (Macmillan,  1922)^  In  the  theory  of  the 
functions  of  a complex  variable.  Other  definitions  of  the  Gamma  func- 
tion are  also  given,  all  of  -which  are  equivalent  to  the  Euler  integral 
in  Eq.  (l)  vhen  the  real  part  of  n is  positive, 
a 

An  integral  with  limits  is  a function  jof  those  limits,  and  not  a 
function  of  the  variable  used  in  the  integrand.  In  Eq.  (2)  the  sub- 
script x in  rx(n)  is  to  be  identified  with  the  upper  limit  of  the  in- 
tegral, not  with  the  x in  the  integrand.  The  student  may  find  it 

helpful  to  revriterEq.  (2)  as  V (n)  = /*  tn-ie-*  dt.  Any  letter  vhat- 

x 0 

ever  may  be  used  in  the  integrand  without  affecting  the  value  of  the 
integral,  provided  its  limits  remain  0 and  x.  Similar  remarks  apply 
to  Eq.  (14) . 
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the  definition  of  the  Gamma  function.  The  curve  is  drawn  for 
n ~ 6.  The  area  under  the  whole  curve  to  infinity  is  the  com- 
plete Gamma  function  T(n),  and  the  area  to  the  finite  abscis- 
sa x is  the  incomplete  Gamma  function  Tx(n). 

not  only  a function  of  n but  also  of  x.  For  consistency 
the  complete  function  r (n)  should  have  the  subscript  ao, 
but  for  brevity  and  by  convention  the  subscript  is  gener- 
ally omitted. 

It  is  instructive  to  plot  the  Integrand  of  Eqs.  (l) 
and  (2)  for  a given  value  of  n,  for  positive  values  of  x. 
The  area  under  the  whole  curve  to  infinity  is  by  definition 
the  complete  function  f(n).  The  area  included  between  0 
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and  the  finite  abscissa  x is  Tx(n).  The  exact  shape  of 
the  curve  varies  with  the  parameter  n,  but  the  relation 
between  the  complete  and  incomplete  functions  is  always  as 
shown  in  Fig.  1.  The  reader  may  agree  that  the  adjective 
incomplete  is  aptly  applied  here;  the  Gamma  function  is 
incomplete  when  the  integration  extends  only  part  way, 
i.e.,  to  some  finite  distance.3 

Other  forms  of  the  above  integrals . --The  above  in- 
tegrals can  take  a variety  of  other  forms  by  change  of 
variable.  In  the  theory  of  errors  an  important  form  is 
obtained  by  setting  x = u2;  then  dx  = 2 u du,  and 

T (n)  = r 2 (n)  = 2 fU  '/X\i*n~1e  “2du.  (3) 

0 

The  student  should  take  careful  note  that  in  the  incom- 
plete Gamma  function  the  subscript  of  T always  refers  back 
to  the  upper  limit2  of  the  integral  in  Eq.  (2);  thus  the 
subscript  x or  u2  in  (2)  is  the  upper  limit  of  the  inte- 
gral in  Eq.  (2)  and  not  of  the  new  form  (3). 

With  n = Eq.  (3)  gives 

r*(-t)  = 2//Xe"x  dx,  (4) 

whence  we  perceive  that  the  normal  probability  integral 
can  be  regarded  as  a table  of  the  incomplete  function 

r*(i). 

Several  other  forms  of  the  integral  for  tne  Gamma 
function  that  are  important  in  application  will  be  found 
in  the  examples  at  the  end  of  this  chapter. 

The  recursion  formula. — The  most  important  proper- 
ty of  the  Gamma  function  is  a recursion  formula.  It  can 
be  found  by  integrating  the  right-hand  side  of  Eq.  (l)  by 
parts,  thus: 


Several  ways  of  evaluating  the  incomplete  Gamma  function  by  series 
and  continued  fractions  will  be  found  in  Exercises  27  and  28. 

4, 

See,  for  example,  E.  B.  Wilson,  Advanced  Calculus,  Ch.  14  (Ginn 
and  Company,  1912) . 


4 


THE  GAMMA  AND  BETA  FUNCTIONS 


r(n)  = /°x n~1  e~x  dx  = (l/n)  /°e-*dxn  (n  > 0) 
o o 

= (l/n)  (U"e-*]“  + f x"e^to} 

= (l/n)  (0  + r(n+  1)  ). 

Accordingly,  there  follows  the  recursion  formula 

T(n+  1)  = n T (n).  (5) 

Thus  F(n+  l)  is  simply  related  to  r(nj.  Eq.  (5)  may  he 
regarded  as  a difference  equation,  and  the  Gamma  function 
as  a solution  thereof. 

By  repeated  applications  of  the  recursion  formula 
it  is  seen  that 

r(n+  l)  = n r(n)  = n r(n-l+l)  = n(n-l)  T(n-l) 

= n (n-  l)  (n-2)  ...  (n-k)  r(n-k).  (6) 

Hence  if  the  Gamma  function  is  tabled  over  a unit  interval 
of  the  argument,  the  function 'for  arguments  outside  this 
range  can  he  found;  for  example 

r(3.37)  = r(2.37+l)  = 2.37 x r(2.37)  = 2.37*  r(r.37+i) 

= 2.37  X 1.37  r(1.37),  (7) 

so  that  if  r(l.37)  is  known,  then!  (3.37)  can  be  evaluated 
hy  a simple  multiplication. 

The  recursion  formula  was  here  derived  from  the 
definition  contained  in  Eq.  (l),  and  it  was  necessary  to 
postulate  that  n > 0;  otherwise  the  term  [xhe-x]0  would 
not  have  dropped  out  in  the  integration  hy  parts.  But  the 
recursion  formula,  it  so  happens,  is  a general  property  of 
the  Gamma  function  however  defined,  and  it  holds  for  nega- 
tive and  complex  values  of  the  argument  as  well  as  for 
positive  real  values. 

The  relation  of  the  Gamma  function  to  factorials. — 

With  n = 1,  Eq.  (l)  gives 
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Hence  if  n is  a positive  integer,  and  if  Eq.  (6)  is  ex- 
tended, the  recursion  formula  leads  to  the  following  sim- 
ple and  important  relation  with  the  factorial: 

r(n+l)  = ii(n-l)  (n-2)  ...  5-4-3-2-1  r(l)  = nl  (8) 

For  values  of  n not  Integral,  r(n+l)  can  he  regarded  as  a 
generalized  nl;  the  practice  is  in  fact  fairly  common. 

Thus  ■fl  is  to  he  identified  with  r("J);  -il  with  rfi");  etc. 

Tables  of  the  complete  Gamma  function. - -Gaus s 5 in 
1815  computed  log  nl  to  twenty  places  and  d In  ni/dn  to 
eighteen  places  in  the  interval  0 Z n ^ 1 hy  steps  of  0.01, 
Degen6  in  1824  computed  log  nl  to  eighteen  decimals  from 
n = 1 to  n = 1200  hy  integral  steps.  Legendre 7 In  1825 
gave  log  r(n)  to  twelve  places  for  1 < n 2 hy  steps 
of  0.001,  together  with  its  first,  second,  and  third  dif- 
ferences. Table  XXXI  in  Tables  for  Statisticians  and 
Biometrlcians,  8 Part  I,  shows  log  r(n)  to  seven  places  In 
the  interval  1 < n < 2 at  intervals  of  0.001,  and  Table 
XLIX  in  the  same  volume  shows  log  nl  to  seven  decimals  for 
every  integer  from  1 to  1000  inclusive.  In  1922  Egon 
Pearson9  published  a table  of  log  r(n)  to  ten  decimals 
from  2 to  1200  inclusive,  by  steps  varying  from  a tenth  to 
a whole  integer.  In  1925  John  Brownlee10  published  a ta- 
ble of  log  r(n)  to  seven  places  by  steps  of  0.01  in  the  . 
interval  1 ? n ^ 50.9.  It  is  now  not  unusual  to  find  ab- 
breviated tables  of  nl  or  its  logarithm  in  tables  for  aids 
in  computation. 


5Carl  Friederich  Gauss,  Werke,  vol.  3,  PP*  161-162.  Gauss’  nota- 
tion is  7iz  for  "Gauss'  infinite  product,"  equivalent  to  zl  He  intro- 
duced ^ z ^ or  the  derivative  In  7iz . 

6Carl  Ferdinand  Degen,  Tabularem  ad  faciliorem  et  breviorem  Probar- 
bilitatis  computationem  utilium  Enneas  (Copenhagen,  1824). 

7Legendre,  Traite  de  fonctions  elliptiques,  vol.  2,  pp.  490-499, 
(Paris,  1825).  A facsimile  of  Legendre’s  tahle  has  been  issued  as 
Tracts  for  Congmters,  No.  4 (Cambridge,  1921). 

8Tables  for  Statisticians  and  Biometricians  comes  in  two  volumes. 
Part  I first  appeared  in  1914,  with  second  and  third  editions  in  1924 
and  1930.  Part  II  appeared  in  1931*  Both  volumes  are  issued  from 
The  Office  of  Biometrika,  University  College,  Gower  Street,  London 

W.C.  1. 

0 

Egon  S.  Pearson,  Tracts  for  Conputers,  No.  8 (Cambridge,  1922). 

10John  Brownlee,  Tracts  for  Computers,  No.  9 (Cambridge,  1923). 
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The  graph  of  the  Gamma  function,. --It  is  easy  to 
evaluate  the  Gamma  function  at  any  positive  integer  with 
Eq.  (8).  For  example,  r (4 ) = 31  = 6 and  r(5)  = 4l  = 24. 
The  curve  is  continuous  on  the  positive  side  of  the  verti- 
cal axis,  as  can  be  proved,11  hence  a graph  like  Fig.  2 is 
not  difficult  to  construct  for  positive  values  of  n.  The 
minimum  between  n = 1 and  n = 2 is  interesting,  and  re- 
quires special  analysis.  We  have  already  noticed  that 


11  See,  for  example,  H.  S.  Carslav,  Fourier’s  Series  and  Integrals 
pp.  132  and  133  in  the  1930  edition  (Macmillan,  1906,  1921,  1930). 
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T (l)  = 1,  and  by  Eq.  (8)  it  follows  that  r(2)  = II  = 1 al- 
so; hence  there  must  be  a proper  minimum  somewhere  between. 
Calculations  have  located  the  position  of  the  minimum  at 
n = 1.46165  21449  68  . ..,  at  which  point  the  Gamma  func- 
tion has  the  value  Q.8856O  51944  ...  . 

For  negative  values  of  n the  integral  in  Eq.  (l) 
is  not  used  to  define  the  Gamma  function;  instead,  other 
definitions  are  used,  such  as  Gauss'  (below)  and  Weier- 
s trass'.  The  recursion  formula  permits  evaluation  of  the 
Gamma  function  for  negative  n in  terms  of  its  values  for  a 
positive  n.  Exercises  5 and  6 at  the  end  give  hints  on 
ways  of  discovering  the  shape  of  the  graph  of  the  Gamma 
function  where  n is  negative.  As  might  be  inferred  from 
Fig.  2,  when  n is  a ne'gative  integer  the  Gamma  function 
suffers  infinite  discontinuities  in  the  nature  of  simple 
poles,  becoming  alternately  positive  and  negative  infinity. 
Accordingly,  the  factorial  of  a negative  integer  such  as 
-II,  -21,  -51,  etc.  is  infinite,  and  its  reciprocal  is 
zero.  This  phenomenon  appears  in  Exercises  15  and  14. 

A very  interesting  chapter  on  the  subject  is  found 
in  Edwards'  Integral  Calculus  (Macmillan,  1922),  Ch.  XXIV 
in  volume  2,  sections  886  and  922  in  particular.  Edwards 
shows  graphs  of 

y = 1/  (n+  1) 


y = 2 ! 2 "/ ( ri  + 1)  (n-t  2) 

y = 3!  3"/ (n+l)(n+2)(n+3) 

y = iuV/(n+l)(n+2)(n+3)(n+4). 


By  continuation,  this  sequence  of  curves  can  be  made  to 
approach  T(n+l)  for  any  n as  closely  as  desired,  where- 
fore the  successive  curves  give  better  and  better  approxi- 
mations to  the  graph  In  Fig.  2 for  finite  values  of  n. 

This  sequence,  continued  indefinitely,  is  in  fact  the  Gauss 
definition  of  r(n+  l)--Gauss ' infinite  product,5  identified 
as  nl 


Tables  of  the  incomplete  Gamma  function. --By  trans- 
formation of  variable,  as  called  for  in  Exercises  52-42, 
several  important  curves  (the  Pearsoh  Types  III  and  V,  and 
others)  can  be  changed  into  the  form  y = x71-1  e~x  , the  in- 
tegrand of  Eq.  (l),  wherefore  the  incomplete  Gamma  function 
constitutes  an  important ' probability  integral.  The  need  of 
tables  of  this  integral  was  recognized  by  Karl  Pearson,  who 
with  the  assistance  of  his  staff  at  the  Biometric  Labora- 
tory, turned  out  in  1921,  after  many  years'  labor,  the 
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Tables  of  the  Incomplete  Gamma  Function.12  A sample  page 
is  shown  in  Fig.  5.  Therein  the  function  l(u,p)  is  tabled 
to  seven  decimals  against  the  arguments  u and  p,  along 
with  second  and  fourth  differences  with  respect  to  both  u 
and  p.  The  function  l(u,p)  is  related  to  u and  p,  and 
through  them  to  the  incomplete  Gamma  function,  by  the  def- 
inition 

l(u,p)  = r„(p+i;/r(p+l)  = /Xx£e~*  dx//°x*>  e~* dx  (9) 

x 0 0 

wherein  u = x//p  + 1 . 

Evidently  for  any  fixed  value  of  p,  l(u,p)  varies  from  0 
to  1 as  x varies  from  0 to  oo;  l(u,p)  is  in  fact  just  the 
fractional  part  of  the  area  under  the  curve  of  Fig.  1 be- 
tween 0 and  the  abscissa  x,  p being  identified  as  n-  1, 
and  u as  x/Vp  + 1 = x//n.  The  advantage  of  tabling 
against  x/Vn  instead  of  against  x directly  arises  from 
the  fact  that  as  n increases  the  curve  of  Fig.  1 flattens 
out,  and  one  is  accordingly  compelled  to  move  to  higher 
and  higher  values  of  x if  he  would  include  always  a speci- 
fied fraction  of  the  area.  The  device  of  using  x//n 
amounts  to  a change  of  scale  that  keeps  the  table  within 
bounds;  when  n is  high,  a moderate  value  of  x//n  corre- 
sp'onds  to  a high  value  of  x.  x//n  is  in  fact  just  the 
upper  limit  x of  the  incomplete  Gamma  function  (Eq.  2)  ex- 
pressed in  units  of  the  standard  deviation  y n of  the 
curve.  (Exercises  1 and  29  will  illuminate  this  point. 

It  may  be  noted  here  that  the  mean  of  the  curve  of  Fig.  1 
lies  at  abscissa  n,  and  the  standard  deviation  of  the 
curve  is  /n.)  The  introductory  material  in  the  Tables  of 
the  Incomplete  Gamma  Function  gives  a history  of  the  work 
and  instructions  for  the  use  of  the  tables,  together  with 
illustrative  examples. 

In  1930  L.  R.  Salvosa  13  published  an  extensive  ta- 
ble of  the  areas  and  ordinates  of  the  Pearson  Type  III 
curve 


12  Tables  of  the  Incomplete  Gamma  Function,  published  in  1922  by  his 
Majesty's  Stationery  Office,  and  reissued  in  193^  by  the  Office  of 
Biometrika,  University  College,  London  W.C.  1. 

13  L.  R.  Salvosa,  Ann.  Math.  Statistics,  vol.  1,  1930;  pp.  191-225. 
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wherein 


y = yo(l  + i’at 

Jo  = ( as)  / e 


~ 1 e-2t/ a 

'/a£  r(4/a2) 


(10) 


the  arguments  being  oc  and  t.  |-ot  is  the  "skewness"  of  the 
curve,  and  t the  abscissa  measured  from  the  mean.  Salvo- 
sa ' s table  can  be  regarded  as  another  form  of  a table  for 
the  incomplete  Gamma  function.  (See  Exercise  46.) 

Tables  of  "chi-square"  can  also  be  used  to  evalu- 
ate the  incomplete  Gamma  function.  Some  of  these  tables 
show 


P(x)  = /"(xs)ik‘VSX  dx8  (11) 

r(ik)22*  x 

tabulated  against  X2  and  k as  arguments.  Other  tables, 
following  Fisher,  show  X2  tabulated  against  P(x)  and  k 
as  arguments.  It  can  readily  be  shown  by  a change  of 
variable  that 


P(X)  = 1 - I\  2(ik)/r(|k).  (12) 

2 X 

(See  Exercises  54  and  42.) 

Another  notation'  for  the  incomplete  Gamma  function, 
used  occasionally  by  continental  writers,  is  (n,x)i  for 
Vx  (n+l).  In  this  symbolism  (n,oo)»  means  nl 

Besides  supplying  tables  for  integrals  under  the 
Type  III  curve  and  other  curves  reducible  thereto,  the  in- 
complete Gamma  function  affords  a method  of  summing  any 
number  of  terms  of  the  Poisson  exponential  limit  without 
approximation,  as  is  demonstrated  in  a later  section  (p.  20). 

Definition  of  the  Beta  function. --The  Beta  function 
will  here  be  defined  by  the  integral14 


14  The  integral  given  here  as  the  definition  of  the  Beta  function  was 

r 1 m-1  n 1-k/n 

studied  by  Euler  mainly  in  the  form  J0  x (1  - x ) dx,  n fixed 
(Art.  345  in  vol.  1 of  his  Institutlones  Calculi  Dlf f erentialis , Pet- 
rograd,  1768),  and  is  sometimes  referred  to  as  "Euler’s  first  integral. 
The  relation  between  the  Gamma  and  Beta  functions,  shown  here  as.Eq. 
(24),  was  given  by  him  in  Art.  27  of  vol.  4. 
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B(m,n)  =/]xri  (1-  x)71-1  dx.  (15) 

0 

B(m,n)  Is  read  "the  Beta  function  of  m and  n. " It  is  tru- 
ly a function  of  m and  n;  that  is,  it  depends  on  both  m 
and  n for  its  value.  With  the  upper  limit  unity,  as  here 
written,  the  Beta  function  is  said  to  be  complete . but,  as 
with  the  Gamma  function,  the  adjective  is  usually  omitted 
except  when  there  is  a possibility  of  confusing  the  com- 
plete with  the  incomplete  Beta  function.  The  incomplete 
Beta  function  is  defined  by2 

B„(m,n)  = x)"--2  dx,  (14) 

0 

the  upper  .limit'  x lying  between  0 and  1.  Evidently  the 
incomplete  Beta  function  B*(m,n)  is  a function  not  only  of 
m and  n,  but  also  of  the  upper  limit  x.  In  the  same  sym- 
bolism the  complete  Beta  function  B(m,n)  should  have  the 
subscript  1,  but  for  brevity  and  by  convention  it  is  gen- 
erally omitted,  just  as  the  subscript  c d is  omitted  in 
writing  r(n)  for  the  complete  Gamma,  function. 

The  student  should  note  that  the  Beta  function, 
complete  or  incomplete  respectively,  has  one  more  argument 
than  the  Gamma  function,  complete  or  incomplete. 

Interchanging  the  arguments  m and  n.--By  change  of 
variable  the  above  integrals  assume  a variety  of  forms4 
and  lead  to  several  interesting  relations,  of  which  per- 
haps the  most  important  is  that 

B(m,n)  = B(n,m)-.  (15) 

This  is  easily  proved  by  replacing  x in  Eq.  (1J>)  by  1 - y. 
If  x = 1 - y,  then  dx  = -dy,  and 

B(m,n)  = (1  - x)*"-*  dx  = - J°(l  - y)m_iy  n~1  dy 

= / V -1  (1-  y^dy  = B(n,m). 

0 

By  the  same  change  of  variable 

Bx(m,n)  = ynJ  (l-  y)m~1  dy 

= J1Vn~1  (1-  y)  ”H‘dy  -f1~XYn~1(l-  y)”~J  dy. 
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That  Is, 

B*(m,n)  = B (m,n)  - B ,_x(n,m).  (l6) 

Or,  in  another  form, 

Bx  (n,m)  = B (ni,n)  - B._x(m,n).  (17) 

Hence  in  the  incomplete  Beta  function  Bx(m,n)  an 
interchange  of  the  two  arguments  alters  the  value  to 
B(m,n)  - That  an  interchange  of  the  arguments 

m and  n has  no  effect  on  the  complete  Beta  function  B(m,n) 
follows  at  once  by  putting  x = 1 in  Eqs.  (l6)  or  (17)  and 
noting  that  B0 (m,n)  = B0(n,m)  = 0, 

These  relations  are  easily  recalled  by  using  a 
plot  of  the  integrand  of  Eqs.  (13)  and  (l4),  such  as  Fig. 

3.  The  effect  of  interchanging  m and  n is  to  produce  a 
mirror  image  of  the  figure,  which  can  be  visualized  as  the 
curve  that  would  be  obtained  by  moving  the  origin  one  unit 
to  the  right,  reversing  the  positive  direction  of  x,  and 
then  looking  at  the  page  through  the  back  side  of  the  pa- 
per. By  visualizing  this  process,  the  relations  exhibited 
by  Eqs.  (l6)  and  (17)  become  evident. 

A recursion  formula  for  the  incomplete  Beta  func- 
tion is  given  in  Exercise  16. 

Unless  m and  n are  equal,  the  curve  y=xr""-z(l-  x )n~1 
will  be  unsymmetrical.  By  varying  the  choice  of  m antf  n, 
a great  variety  of  forms  can  be  obtained.  The  student 
will  profit  by  carrying  out  the  calculations  and  plotting 
the  curves  suggested  in  Exercise  2. 

Two  other  useful  forms  of  the  Beta  function. --Let 
x = z2/(l+z2),  1 - x = 1/(1+ z2),  dx  = 2z  dz/(l+z2)2.  By 
introducing  this  change  of  variable  into  Eq.  (13 ) we  find 
that 

B (m,n)  = 2 /“ z2*  ~ 1 {l  + z2)"n,-ndz  = B(n,m). 
o 

In  particular,  if  m is  replaced  by  £ and  n by  i[n 
this  relation  gives 

B(i[n  - 1] , i)  = 2 / ” (1+  z s)~in  dz  - /“  (1  + z2)-+dz,  (19) 

0 “°° 

which  will  be  found  handy  for  certain  kinds  of  calcula- 
tions. 
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grand  in  the  definition  of  the  Beta  function.  The  curve  is 
drawn  for  m = 5,  n = 3*  The  area  under  the  whole  curve  from  0 
to  1 is  the  complete  Beta  function  B(m,n),  and  the  area  Between 
0 and  x < 1 is  the  incomplete  Beta  function  B(m,n). 

Another  useful  form  is  obtained  by  the  substitu- 
tion x = sin2 9,  dx  =.  2 sin  9 cos  9 d 9,  whereupon  the  orig- 
inal form  shown  in  ' Eq.  (13)  goes  into 

B'(m,n)  = 2 sin2”1  “-*9  cos  2n~1y  d9.  (20) 

o 

In  particular, 

= 2f2n  a?  = k. 

a relation  that  vill  be  needed  shortly  (see  Eq.  26). 


(21) 
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The  relation  between  the  Gamma  and  Beta  functions.-" 
From  the  second  form  of  the  Gamma  function,  Eq.  (3), 

r(m)r(n)  = 4/ <r  u2m-ie-t/2du-/°C  v2n“ie”,y2dv 
0 0 

= 4/  / u.2m-lv2n-xe~  (u2+v2)du  dv.  (22) 

o o 

It  is  permissible  to  change  the  product  of  two  Integrals 
to  a double  integral,  and  vice  versa,  when  neither  the 
limits  nor  the  integrands  of  either  integral  depend  on  the 
variable  in  the  other.15 

Now  consider  the  surface 

Z = U2m-lv2n-le-(u2+v2 ) (23) 

erected  on  the  uv  plane.  u2m~1v2n~1e~ (w2+v 2)du  dv  can  be 
considered  the  volume  of  an  elementary  rectangular  prism 
erected  on  the  base  du  dv  and  having  altitude  z,  as  shown 
in  Fig.  5.  The  integral  written  in  Eq.  (22)  calls  for.  the 
summation  of  the  volumes  of  all  these  elementary  prisms, 
by  which  process  one  finds  the  total  volume  under  the  sur- 
face generated  by  Eq.  (23).  The  summation  (integration) 
can  be  carried  out  by  introducing  polar  coordinates  in  the 
uv  plane . Let 

u = r cos  qp,  v = r sin  qp. 


Then  the  element  of  area  du  dv  in  Eq.  (22)  is  to  replaced 
by  r dr  dqp.  A quadrant  of  the  surface  described  by  Eq. 
(23)  can  be  mapped  out  by  allowing  qp  to  vary  from  0 to  j-rc 
and  r from  0 to  00,  That  the  integral  in  Eq.  (22)  actually 
exists  (i.e.,  approaches  a definite  value  as  the  upper 
limits  for  u and  v are  made  larger  and  larger)  can  be  dem- 
onstrated by  showing  that  the  volume  under  the  surface  of 
Eq.  (23)  outside  the  radius  R can  be  reduced  to  any  de- 
sired value  however  small  simply  by  making  R big  enough. 
Eq.  (22)  then  leads  to  the  following  equalities: 


f(m)  f(n)  = 4/  V2m-lv2n-le-rzv  dr/2  sin2*1"1  qp  cos  2n~1  <p  d qp 


= 2 /°°  r2(n+»)-le-r2dr  # 2 j 2n  s±n2n-l(p  cos  2m-l(p  d 


= T(nfm)B(n,m) . 


15E.  B.  Wilson,  Advanced  Calculus, Art . 143  (Ginn  and  Company,  1912). 
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Fig.  5»  The  elementary  prism  of  height  z and  base  du  dv,  cut 
from  the  surface  z = u v ^"^e  + v'  . The  volume  under  the 
whole  surface  is  the  sum  of  the  volumes  of  all  such  prisms. 


by  Eqs,  ( 5 ) and  (20).  Thus  we  have  discovered  that 


B (m,n) 


Lfelrjni 

r (m+n) 


(24) 


Another  proof  of  this  relation  -will  be  found  in  Exercise 
11  at  the  end. 

To  find  T(?) . — At  the  close  of  the  paragraph  on 
other  useful  forms  of  the  Beta  function  it  was  found  that 

= *.  (See  Eq.  21,  p.  1?) 


Now  with  m = n = § In  Eq;  (2k)  we  have 


r(j)r(j) 
r(D  * 


and  since  r(l)  = 1 It  follows  that 


(25) 


r(i)=v/it.  (26) 

Eq.  (3)  then  gives 

r(i)  = 2fWs-y2&j  = f°“ e~yS  dy  = /it  . (27) 

0 —go 

The  same  result  is  derived  in  Exercise  25.  By  means  of  the 
integral  just  written,  the  student  should  satisfy  himself 
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that  the  normal  curve  when  written  in  the  form 

y 7=r  e~*s/2a2  . (28) 

q/27t 

is  normalized,  I.e.,  that  the  area  under  the  curve  of  y 
plotted  against  x is  unity,  and  that  it  would  not  be  so  if 
the  constant  factor  were  other  than  l/a^/Sht. 

Tables  of  the  incomplete  Beta  function. --The  ratio 
Bx (m,n)/B (m,n)  is  the  fractional  part  of  the  total  area 
under  the  curve  y = xm~l  (l-  x)n~l  (Fig.  J>)  between  0 and 
the  abscissa  x < 1;  and  1 - B X(m,n)/B(m,n)  or  B i_x  (n,m) 
/B(m,n)  is  the  fractional  part  remaining  between  x and  1, 
This  curve,  which  is  the  graph  of  the  integrand  of  the 
Beta  function,  constitutes  another  important  frequency 
curve,  the  Pearson  Type  I,  and  by  change  of  variable,  the 
Pearson  Types  II,  VI,  and  VII  also;  hence  the  Incomplete 
Beta  function  constitutes  another  important  probability 
integral.  The  immense  task  of  providing  such  tables  was 
finally  realized  in  19354  with  the  publication  by  Karl 
Pearson  and  his  staff  of  the  Tables  of  the  Incomplete  Beta 
Function  16  --an  even. greater  undertaking  than  the  Tables 
of  the  Incomplete  Gamma  Function.  A sample  page  is  shown 
in  Fig.  6.  The  function  tabled  is 

= B*(p,q)/B(p,q),  (29) 

the  arguments  being  p,  q,  and  x.  The  introduction  gives  a 
history  of  the  incomplete  Beta  function  and  illustrations 
of  its  use.  The  complete  Beta  function  B(p,q)  is  listed 
in  the  headings  of  the  columns,  so  that  if  desired, 

B*(p,q)  can  be  found  by  multiplying  Ix(p,q)  by  B(p,q). 

Besides  supplying  tables  for  integrals  under  cer- 
tain frequency  curves,  the  Incomplete  Beta  function  sums 
terms  of  the  point  binomial  without  approximation,  as  is 
demonstrated  in  a later  section  (p.  18) . 


16TableB  of  the  Incomplete  Beta  Function,  issued  "by  the  Office  of 
Biometrika,  University  College,  Gower  Street,  London  W.C.  1.  Mention 
should  he  made  of  Tracts  for  Computers  No.  7y  hy  H.  E.  Soper,  entitled 
The  Numerical  Evaluation  of  the  Incomplete  Beta  Function,  which  will 
he  found  useful  in  the  absence  of,  or  outside  the  range  of,  the  Tables 
of  the  Incomplete  Beta  Function. 
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Fig.  6.  A sample  page  in  the  Tables  of  the  Incomplete  Beta  Function. 
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To  sum  any  number  of  terms  of  the  point  binomial 
without  approximation. --Let  S denote  the  sum  of  the  terms 
of  the  point  binomial  (q  + p)n  as  far  as  the  term  in  p t, 
which  is  to  define 


S = qn  + nq n~*  p + (£) qn”2  p2  + . . . -S-  (^Jq”"*  pt 


= 2 e)qn-rpr.=  q -r  pr 

r—0 


7-n  'r'  ‘ * (n-r).’r! 


(50) 


q and  p are  to  be  complementary;  i.e.,  q = 1 - p,  whence 
dq/dp  = -1.  Then  if  S be  differentiated  with  respect  to 
p,  the  result  is 


ni 


dS  *_  t 

dp  ~ “ 2q  (n-r-l)  irl 


q n-r-1 pr  + £ 


ni 


i-l 


(n-r) I (r-l) 1 


q n-r^r-l 


n • — a n-t-1  t 

, q P . 


(n-t-1) It  I 


. J Hi n-s-ls 

+ sTo  (n-s-l)ls!  q p 

(51) 


This  simple  result  is  produced  by  the  cancellation  of  plus 
and  minus  terms  in  the  two  summations;  for  every  term  in 
the  first  sum,  except  for  r = t,  there  is  a term  of  oppo- 
site sign  in  the  second  sum. 


Now  when  p = 0,  S = 1,  whence  from  the  last  equa- 

tt  *-*-* 


tion, 
which  gives 


dx 


(32) 


S = 1 


n! 


(n-t-1) It!  o 


fyxt  (1  - x ) n~t~1  dx.  (55) 


Moreover,  when  p = 1,  q = 0 and  S = 0 if  t < n,  hence  the 
factorials  must  satisfy  the  relation 

A*(l d,  . . 

0 H • 


(3^) 
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This  relation  is  equivalent  to  Eq.  (24),  hut  the  evaluation 
of  the  factorials  in  terms  of  the  Beta  function  comes  in 
here  as  a by-product.  It  follows  at  once  that 


/V(  1 - x)  dx 
S = 1 - -2— 

//x‘(l  - x)"-*-1  dx 


(35) 


Thus  S turns  out  to  be  unity  diminished  by  the  ratio  of  an 
incomplete  to  a complete  Beta  function,  and  no  approxima- 
tion has  been  introduced.  Unfortunately,  the  use  of  the 
incomplete  Beta  function,  even  with  tables,  is  often  labo- 
rious. Indeed,  Eq.  (35)  can  and  has  often  been  used  as  a 
means  of  evaluating  the  integral  /Q  x*  (l-  x)n-t~1  dx  by 
summing  terms  of  the  binomial;  such  was  in  fact  the  use 
that  Bayes  made  of  it  in  1763  (Trans.  Royal  Society  for 
the  year  1763,  p.  396).  Laplace  later  derived  the  same 
result  with  the  opposite  purpose,  viz.,  to  sum  terms  of 
the  binomial  (Theorie  analytique,  1823;  p.  151 )• 


The  derivation  followed  here  was  published  by  Laz- 
arus 17 ; it  is  simpler  than  any  of  the  others  that  I have 
seen. 


To  sum  any  number  of  terms  of  the  Poisson  exponen- 
tial limit  without  approximation. --Let  S denote  the  sum  of 
the  terms  of  the  Poisson  exponential  limit  as  far  as  the 
term  in  mL  Then 


S = e-ro{  1 + m + m2/2I  + ...  + m* /t I } . (36) 


Differentiation  with  respect  to  m gives 

H = -S  + e ~m{  1 + 2m/2!  + 3m2/3 1 + ...  + t m'1  /tl  } 

= -e~nmt/t\  (37) 


17This  method  of  deriving  Eq.  (35)  was  first  published  by  William 
Lazarus  of  Hamburg  in  the  Journal  of  the  Institute  of  Actuaries,  vol. 
15,  1870;  pp.  245-257;  page  251  in  particular.  Lazarus  credits  this 
derivation  to  one  Dr.  Landi  of  Trieste.  The  result  contained  in  Eq. 
(35)  was  first  found  by  Bayes;  Trans.  Royal  Soc.  for  the  year  1763, 
page  396.  This  paper  is  available  in  a booklet  entitled  Facsirniles 
of  Two  Papers  by  Bayes,  with  commentaries  by  E.  C.  Molina  and  W.  Ed- 
wards Deming,  published  by  the  Graduate  School  of  the  Department  of 
Agriculture,  1940, 
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all  the  terms  in  the  braces,  except  for  the  last,  being 
cancelled  by  opposite  terms  in  S.  Now  when  m = 0,  S = 1, 
so 

S w 

f as  = (l/tl)  / x*  e~x  dx 
1 0 

m 

S = 1 - (l/t I ) / x*e~*  dx.  (58) 

0 

Moreoever,  when  m = 00 , S = 0,  and  it  follows  that 

00 

f xte~x  dx  = t!  (39) 

0 

This  relation  is  equivalent  to  Eq.  (8),  but  the  identifi- 
cation of  the  complete  Gamma  function  with  the  factorial 
comes  in  here  as  a by-product.  It  follows  at  once  that 

m x 

S = 1 - / x*e“*  dx/  / x*  e“x  dx.  (40) 

0 0 1 

Thus,  the  incomplete  sum  of  the  Poisson  exponential  limit 
is  expressible  in  terms  of  an  incomplete  Gamma  function. 
Turned  around,  the  incomplete  Gamma  function  is  expressi- 
ble as  a sum  of  terms  of  the  Poisson  limit.  Eq.  (40)  was 
first  published  by  E.  C.  Molina  in  1915  (Amer.  Math.  Month- 
ly, vol.  22,  p.  223,  1915).  It  is  usually  derived  by  in- 
tegration by  parts,  but  Lazarus'  method  seems  simpler. 


EXERCISES 


1.  Plot  the  curve  y = x 71-1  e~x  for  n = 1,  2,  5,  6, 
10  between  x = 0 and  x = 8.  The  areas  under  these  curves 
to  infinity  are  respectively  T(l)  = 1,  T(2)  = 1,  P(5)  = 2, 
r(4)  = 51  = 6.  Notice  that  as  x ->a>  the  slopes  of  all  the 
curves  become  zero.  If  n = 1 the  curve  starts  from  the 
point  (0,1)  with  slope  -1;  if  n = 2 the  curve  starts  from 
the  origin  with  slope  +1.  If  n > 2 the  curve  starts  from 
the  origin  with  a horizontal  tangent,  making  contact  of 
order  n - 2.  A.t  infinity  all  the  curves  make  contact  of 
infinitely  high  order  with  the  x axis. 

2.  (a)  Plot  the  curve  y = x7”-^  (l-  x)”_ 1 for  m, 

n = 1,  2;  2,  1;  2,  3;  5,  2;  J,  2;  2,  \ between  x = 0 and 
x = 1. 
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(b)  Determine  the  numerical  values  of  the  areas 
under  these  curves,  and  check  roughly  with  the  graphs.  No- 
tice that  when  n = m the  curves  are  symmetrical,  and  when 
n ^ m the  curves  are  lopsided;  if  n < m the  bulk  of  the 
area  is  thrown  to  the  right,  and  if  n > m the  bulk  of  the 
area  is  thrown  to  the  left.  The  order  of  contact  with  the 
x axis  is  m - 2 at  x = 0 • and  n-2  at  x = 1. 


3.  Given  r(~)  = y/n;  find  by  the  recursion  formula 
the  values  of  f(|),  f(|) , f(|) , f(|) . Show  that  (n  + -)  l 


= f(n  + f ) = 


ni  — 


2 n+1 


(2n  -l-  I) 


cc 

4.  By  the  integral  r(n)  = / x”~-ze-xdx,  show  that 

0 


r(n)  increases  without  limit  as  n + 0 from  the  positive 
side.  Proof:  Near  the  origin,  e~x  can  be  replaced  by  1, 


so  that  the  integrand  becomes  x n~l . The  integral  /°xn  1 dx, 

0 

0 < a,  is  finite  as  long  as  0 < n;  but  no  matter  how  small 
a is,  this  integral  can  be  made  greater  than  any  preas- 
signed number  however  large  by  making  n small  enough.  In 
other  words  the  integral  approaches  +<o  as  n ->■  0 from  the 
positive  side.  Hence 


Lim  r(n)  = +<» 

n -+  + 0 

(See  Fig.  2.) 

5.  (a)  By  the  recursion  formula  show  that  r (--jp)  = 
-2r(|)  = -2/ti  . 

(b)  r(-n  + —)  = /Jt(-2)  Vl-3.5.  . .(2n-  1). 

Notice  that  r(-n  + alternates  in  sign  and  decreases  in 

magnitude  as  n increases.  This  formula  evaluates  the  Gam- 
ma function  at  the  negative  half-integers;  see  Fig.  2. 

6.  Show  by  the  recursion  formula  that  r(-^)  = 

4r(f);  r(-i)  = -8r (■£),,  f(-i)  = -52r(g),  ....  r (-f ) = 
,T.(j)r(|),  r(-|)  = -(|)r(|),  r(-||)  = ...  . so 

p(n)  approaches  -«  as  n approaches  0 from  the  negative  side 
or  -1  from  the  positive  side. 

7.  Show  that  if  ni  be  identified  with  F(n+l),  then 

01  = 1. 
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8.  Show  that 


oo  , 2 2 

bx-a  x dx 


/ e 


b2/  4 a2 


2 " Tb2/4a2  b = 0 

00  2 2 e 

/ e6l-a  * dx  = — 2|a|  [»*  + r62/4a2  (I)],  b > 0 


By  addition 

2 / e-a^*2cosh  bx  dx  = / e bx-a2x 2 dx  = — 


b2 /4a2 


b = 0. 

o o |a-i  > 

Given:  cosh  x = i(e*  + e-*). 

9.  /^t”--* e~ 4 ”dt  = f^t271-!  e~tn  dt  = — . 
oo  n 

/ tw”-ze“tn  dt  = (— )T(— ),  which  becomes  n2/n  if  n = 2m. 
o n n 

(Kramp,  Analyse  des  refractions,  Strasbourg,  1799) • 


10.  r( 


n-1 


n-3  n-5  n-7 

2*2*2 


n-3  . n-5  n-7 

2*2*2 


1 . 5 

2 2 


/ 71 


/n 


(n-2) 


= [*(n-3)].' 


6 4 

2 ’ 2 


n even 


n odd 


2 2 
n even 

2 

2 

n odd 


By  means  of  these  verify  your  answers  for  Ex.  J>. 

I I . By  trigonometric  substitution,  show  that  the 

Euler  integral  ( 1-x) n~1  dx  defining  B(m,n)  is  equiv- 

in  o 

alent  to  2 f2  3ln2m~1  0 cos2n_J  0 d 0.  Integrate  this  by 
o 

Wallis’  formula  and  show  that 

(m-l) ! (n-l) ! . r(m)r(n) 

B(m,n)  = B(n,m)  = (nH.n.i)—  - rf^n)-  ' 


which  is  Eq.  (24). 
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1 

r(in) 

r n— 1 lx 

‘ 2 ' 2j 

’ r(S§i)  r(|) 

1 

n-2 

n-4 

n-6 

6-4-2 

2n-2(S|2;)2 

n 

n-3 

n-5 

n-7 

5-3-1 

?t(n-2) 

1 

n-2 

n-4 

n-6 

...  5-3-1 

(n-2).' 

2 

n-3 

n-5 

n-7 

4*2 

2n-2(Sz3;)2 

n even 


n odd 


13.  The  number  of  possible  combinations  of  n arti- 
cles taken  r at  a time  is  nl /(n-r ) ! r J , vhich  is  often  ab- 
breviated (5),  as  in  Eq.  (30), 

(a)  Show  that  if  the  factorials  in  this  expres- 
sion be  identified  with  Gamma  functions,  then  (?)  vanishes 
when  r is  any  negative  integer  or  any  positive  integer 
greater  than  n (a  fact  already  mentioned  in  the  text). 

(b)  Show  that  .(g)  = (§')  = 1. 

(c)  Show  that  2r(?)prq"  = np  where  the  sum 

is  to  be  taken  over  all  integral  values  of  r. 

n 2 

I".  Show  that  ri0  B(n-r+l,  P+1)  = 2”(n+l)’  or  that 

(g)  + (£)  + (“)  + •••  + (“)  = 2".  Hint:  Expand  (j?  + £)" 

by  the  binomial  theorem  and  note  that  the  sum  of  the  terms 
is  unity.  The  summation  could  as  well  be  taken  over  all 
integral  values  of  r,  positive  and  negative;  see  the  pre- 
ceding exercise. 

15.  Show  that  B(m, n)  = f1  - ^ dx.  Hint:  re- 

0 (1+x)*+n 

place  x in  B(m,n)  = 1-x)  n“-*  dx  by  z/(r+z)  and  then 

0 

by  l/(l+z)  and  get  B(m,n)  = / z m~1(l+z)  ~m~n  dz  and  B(m,n) 

o 

= / .-i(l+z)—«  dz  . Add.  these  and  obtain  B(m,  n)  = 


1*7,™  1 + zn-]  , 1 r1zm~1  + zn~1  1 spz m~1  + z n~l 

-j  — dz  = —j dz  + -f 

° (l+z)"+n  2°  (l+z)m+n  21  (l+z)*+n 


dz 
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= a/ 


1 A z"1-1  + z 


2‘°  (l+z)”,+n 


dz  + -/ 

O 


1 A ' v ” ~ 1 + v™-1  (z  = l/v  in  the 


dv 


(1+v) 


m+n 


last  integral 


= / 


xm-l  + xn-l 


dx . 


(1+x) 


i + n 


(6.  Prove  the  following  recursion  formula  for  the 
incomplete  Beta  function-- 

I*(p,q)  = xlx(p-l,q)  + .(1-x)  Ix(p,q-1) 

Here,  as  in  the  Tables  of  the  Incomplete  Beta  Function, 

I*(p,q)  = B (p,q)/B(p,q)  = - -y  fXxi>~1  ( 1-x)  q~1  dx. 

* * £Hp,q;  o 


Solution  due  to  Mr.  Martin  Katzin, 

Naval  Research  Laboratory,  April  1942. 

The  problem  is  to  prove  the  recursion  formula 


I_(p,q)  = xl  (p-l,q)  + (l-x)l  (p,q-l) 


(1) 


Start  with 


- Ifcf  ■ f dt  (2) 


then  write 


, . _ r(p)r(q)  _ (p-pr(p-i)r(q) 

B P'q  " r(p+q)  (p+q-i) r(P+q-i) 


" ^TB(p'q-v 


(3) 


whence 


(p+q-l)  B(p,q)  = (p-1)  B(p-l,q)  = (q-l)  B(p,q-l)  (4) 
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Also 

B~ ( p - 1 , q ) = fXtp~2  ( l-t)q~ 1 dt 


t P-1 


—(l-t)f-1  | +^jfxtp~1  d-t) 

1 | p-l  o 


q-2 


dt 


±t>-l  ( 1 -x)q~1  q-1 

p-l  + p 


+ Bx  (p, q-1 ) 


(5) 


whence 

(p-l)  B^  (p-l, q)  = (q-1)  Bx  (p, q-1 ) + xp-1  ( 1-x )q~1  (6) 

Divide  these  terms  by  the  three  members  of  Eq.  (4)  and  get 


Ir(p-l,q)  - I*(p,q-1)  = 


X*-1  (l-x)^-1 

(p+q-i )B (p, q ) 


(7) 


Now  work  with  Bx(p,q).  We  need  something  different 
from  any  reduction  yet  accomplished.  Find  now  that 


B (p,q)  = fX  tP  ^ (l-t)*7 

o 


dt 


t* 


(l-t)*-1 


q-if* 

P 0 


+ -^fX  tp(l-t)q~1  dt  (8) 


whence 


■ = Bx(p,q)  - ^fX  dt 


(9) 


It  is  of  course  true  that 


(q-1 )fX  tp(l-t)q  2 dt  = + (q-1 )fX  -t  • t*  J(l-t) 


q-2 


-1 


q-2 


dt 


= (q-l)  /*  (1-t-l)  dt 


.*  .P-l 


= (q-1)  f*  ±{l-t)q  1 dt  - (q-1)  r t^‘" (l-t) 


x t-1 


q-2 


dt 


® (q-l)  Bx(p,q)  - (q-1)  Bx(p,q-l) 


(10) 


So 


p(l-x)q  1 = P B (p*-q)  + (q-1)  B (p, q)  - (q-l)  B (p,q-l)  (ll) 
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or 

xft(l-x)  q~1  = Br(p,q)  _ q-1  Bg(p,q-l) 

(p+q-l)  B(p,q)  B ( P , q ) P+q-1  £(p,q) 

= I*(p,q)  - Ix(p,q-l)  (by  Eq.  (12 

Also,  from  Eq.  (7)  we  have 

and  .by  comparison  it  is  seen  that 

I*(p,q)  = x Ix(p-l,q)  + (l-x)  Ixfp,q-l),  Q.E.D. 

17.  In  the  analysis  of  variance  occurs  the  fre- 
quency curve 


y dv  = 


% k 1 

2 k’2  k 


1.  b" 

ft  2 « 


B(ik’,  ±k")  (k"  + k-w2)2 


2 \y(  k'  + k 


7) 


dw 


for  the  distribution  of  w,  which  stands  for  the  ratio  of 
two  estimates  of  cr.  By  transforming  to  x,  show  that  if 
P(v)  denotes  the  fractional  part  of  the  area  lying  beyond 
a given  abscissa  w,  then 


P(w)  = 1 - 


fx  xP_i (l-x) Q-1  dx 


B(p,q)  0 

- 1 - ix(p,q)  = Ii_x(^p) 

wherein  x stands  for  k 'w2/(k"+k 1 w2 ) , and  p = Ik',  q = |kn 
Thus  the  evaluation  of  P(w)  can  be  made  to  depend  on  an 
incomplete  Beta  function. 


18.  The  Type  III  curve  Is  sometimes  written 
y = y (1  + x/a)^  e~*>X//a  in  order  to  put  the  zero  of  x at 

the  mode  (maximum) . Select  a point  B lying  at  the  dis- 
tance b from  the  finite  end  of  the  curve,  and  show  that 
the  fractional  part  of  the  area  lying  beyond  the  point  B 
is 

/“  V*  e~v  dv//“  V*  e~v  dv  = 1 - l(u,p) 
bp / a f 0 
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if  u = bp/a(l+p)2  in  the  notation  of  the  Tables  of  the  In- 
complete Garmna  Fane t ion  ( p . 9 ) • 

19.  Prove  that  ftxm~1(a-x)  n~1  dx  = &m+n~1  B f / (m,n), 

a > 0.  - 0 

20.  The  volume  bounded  by  the  xy  plane,  the  plane 
x+y  = 1,  and  the  surface  z = x ym~1J  Is  T ( /)  T(m)/ 

r(/  + m + 1)  . 

21  . (a)  The  integral  fff  x^'1  jm~ 1 z n~1  dx  dy  dz 
taken  over  all  positive  values  of  x,  y,  z such  that 
x+y+z  ^ 1 is  equal  to  T(/)  Tfm)  r(n)/lY/  + m + n + l). 

(b)  The  integral  fff  x^1  jm~ 1 zn~1  dx  dy  dz 
taken  over  all  positive  values  of  x,  y,  z such  that  x/a 
+ y/b  + z/c  ^ h,  is  equal  to  a / b m c n h^+m+ri  V{f)  r(m) 

r(n)/It/  + m + n .+  l).  Hint:  Transform  the  integral  of 

part  (a)  by  putting  x,  y,  z = ah£,  bhrj,  chC. 

(c)  The  integral  fff  x^"1  jm~1  zn~1  dx  dy  dz 
taken  over  all  values  of  x,  y,  z lying  -within  the  positive 

octant  of  the  surface  (x/a)^  + (y/b)*7  + (z/c)r  ^ h is 
equal  to 

ht/t  + m/q  + n/r  ^b'c'  VjJ/v)  r(m/q)  r(n/r) 

pqr  f(//p  + m/q  + n/r  + l) 


Hint:  Transform  to  an  Integral  like  that  in  part  (a)  by 

putting  £h,  7)h,  Ch  = (x/a/>  , (y/b)^  , (z/c)r  . 


(d)  By  the  results  of  part  (c)  show  that  the 

volume  of  the  ellipsoid  (x/a)2  + (y/b)2  + (z/c)2  = 1 is 

4 43 

— 7i  abc  and  hence  that  the  volume  of  a sphere  is  — n r . 

[Hint:  Put  h = 1;  p = q = r = 2;  / = m = n = 1 . ] 


(e)  By  a similar  artifice,  show  tnat  the 
volume  of  the  hypocycloid  (x/a Y^°  + (y/b)  2/3+  (z/c)  2^3 
= 1 is  4it  abc/55- 


(f )  The  center  of  gravity  of  the  positive  oc- 
tant of  the  ellipsoid  of  part  (d)  lies  at  x,  y,  z = Ja/8, 
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3b/8,  3c/8.  Hint:  In  part  (c)  put  / = 2,  i = n = 1, 

h = 1,  p=q=r=2to  find  the  x coordinate  of  the  cen- 
ter of  gravity. 

22.  (a)  The  n-fold  integral  ff  . . . f x™1'1  x 2m2~1.  . . 
x^n'1  dxj^dxs.  . .dxn  taken  over  all  positive  values  of  x1} 
x2,  . . . , x n such  that  x\  + x2  + ...  x n + 1 1,  is  equal  to 
rfajJrCma) . . . rfm^J/rfm.,.  + m2  + ...  + mn  + 1).  [This  is 
an  easy  extension  of  part  (a)  of  the  preceding  exercise.] 


(^/aj  + 
niu 

a"  r(|n+l) 
xf  + x|  + 


(h)  The  volume  of  the  n-dimensional  ellipsoid 
2 2 

(x2/aa)  + ...  + (xn/an)  = 1 is  V = axa2  ... 

, and  the  volume  of  the  n-dimensional  sphere 

2 2 t^V1* 

..  + x~  - r is  — — . Hint:  The  volume  V 

” r(|n+l) 


of  the  n-dimensional  ellipsoid  is  by  definition  2n  times 
the  integral  of  axidx2  ...  dxn  taken  over  all  positive 
values  of  xx,  x2,  ...,  xn  such  that  (xi/aj.)2  + (xg/a2)2 
+ ...  + (xn/an)2  1 1.  By  the  change  of  variable  (Xj/at  ) 
= ^ and  dx^  = Ja^“2  d^j,  the  volume  V is  seen  to  be  2n 

£Lj^8»  ^ • • • ^*7i 

times  the  integral  of  y-  d£id£2  ...  d£n 

2"(«iS«  •••  ?n)2 

taken  over  all  positive  values  of  £2,  . ..,  £n  such 

that  + Z2  + ...  + Zn  = 1.  The  result  of  part  (a)  with 
m1=m2  = ...  =mn  = -J  applies  at  once,  giving  V = aj_a2 
in 

7l2 

...  a_  — as  written  above. 

r(±n+l) 


23  . (a)  The  volume  of  the  thin  spherical  shell  ly- 
ing between  the  pair  of  n-'"M  1 ^ 2 j_  ^ 2 

in 


•dimensional  spheres  xf  + xg  + 


+ x = (r  + |dr)  is 


2n‘ 


r(in) 


rn-idr.  [Hint:  simply  differen- 


tiate the  volume  of  the  sphere  just  found.]  This  result  was 
used  by  Helmert  in  developing  the  distribution  of  the 
standard  deviations  in  samples  of  n drawn  from  a normal 
population.  See  Czuber's  Beobachtungsf ehler  (Teubner, 

1891)  on  pages  147-150.  The  student  should  observe  that 
this  result  reduces  to  4rcr2  dr  for  the  volume  of  a 
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spherical  shell  in  three  dimensions,  and  to  2nr  dr  for  the 
area  of  a circular  ring  in  two  dimensions. 


(b)  The  volume  of  the  thin  shell  lying  between 
the  pair  of  n-dimensional  ellipsoids  (xi/ai)2  + (x2/a2)2 


+ (xn/an)2=  (s  + Jds]P  is  aia2 


r(in) 


Sn‘%3. 


This  result  is  useful  in  developing  the  distribu- 
tion of  the  standard  deviations  in  samples  of  n drawn  from 
a normal  population.  See  page  126  in  Deming  and  Birge ’ s 
Statistical  Theory  of  Errors  (The  Graduate  School,  Depart- 
ment of  Agriculture). 


24.  The  area  of  a sphere  of  radius  r is 
{2Tc2ll/r(|-n)  }vn~1 . Hint:  simply  take  off  the  dr  from  the 

result  of  Exercise  23.  The  student  should  observe  that  in 
three  dimensions  this  result  reduces  to  4rcr2  for  the  area 
of  a sphere,  and  in  two  dimensions  to  2 7irfor  the  circum- 
ference of  a circle. 


25.  Euler  in  1768  showed  that 
oo  yn~l  . 

f 7— — dy  = it/sin  nir  if  0 < n < 1. 

0 1 + y 

Use  this  result,  and  put  x = y/(l+y)  in  Eq.  (13)  to  prove 
that 

(a)  T(n)r(l-n)  = B(n,l-n)  = 71/sin  nit 

if  0 < n < 1. 


(b)  Set  n = \ and  prove  that  r(-|)  = /71,  whence 

fCCex  dx  = -J  /it. 

0 

(c)  Prove  that  r(J-)  = n/T(i) 

r(f)  = /2  K/r (i) 
r(f)  = 2*/r(j) 

Etc. 

(d)  Prove  that  r (l/n)r (2/n)r (3/n)  ...  r(n-l/n) 

= (2n)i(»-D/  /n. 
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Hint:  Set  A equal  to  the  left-hand  side.  Then 

A2  = f (l/n) T(n-l/n)  . r(2/n)r(n-2/n)  r(n-l/n)r(l/n) 

= {7t/sin  7t/n}{n/sin  27i/n}  ...  {7i/sin  n-1  7t/n} 

by  part  (a).  Then  let  0 -►  0 in  the  following  trigonomet- 
ric identity  (Hobson:  Plane  Trigonometry) 

sin  nfl/sin  6 = 2n“^sin(  0 + 7t/n)sln  (0+27r/n)  ...  sin(6  + n-l  n/n) 

with  the  result  that  A2  = (27i)  n~1/n,  which  is  equivalent  to 
the  theorem  stated. 

26-  Show  that  the  relation  between  the  Gamma  and 
the  Beta  functions  (Eq.  24)  can  be  found  without  recourse 
to  integration  in  polar  coordinates. 

Hint:  (Called  to  my  attention  by  Mr.  Arnold 

Frank)  Evaluate 

I = /°°  f e~x(y+l  )xm+n”-*y  dx  dy, 
o o 

first  with  respect  to  x,  and  then  with  respect  to  y,  to 
find 

I = T(m  + n)B(m,n)  (Let  z = x(y  + 1)}. 

Now  evaluate  the  same  integral,  first  with  respect  to  y, 
and  then  with  respect  to  x to  find 

I = f(m)r(n)  {Let  z = xy } . 

The  order  of  integration  is  immaterial,  whereupon 
these  two  results  can  be  set  equal,  and  the  relation  be- 
tween the  Gamma  and  the  Beta  functions  follows.  That  the 
order  of  integration  is  immaterial  can  be  justified  by  a 
theorem  due  to  de  la  Vallee  Poussin.  This  theorem  states 
that  the  values  of  certain  infinite  multiple  integrals  are 
independent  of  the  order  of  integration. 

27.  Expansion  of  the  incomplete  Gamma  function  in 


series . 
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(a)  Expand  e~x  in  powers  of  x and  get 

T (n)  = fxxn~1 e“xdx 
X 0 

= {xn/n}{!  - - — — x + n 


n 4-  1 


n + 2 


n + 3 


x3  + . 


This  series  is  convergent  for  all  values  of  x,  hut  is  not 
convenient  for  numerical  calculation  unless  x is  small. 

(b)  By  integrating  by  parts  show  that 

fXxn~^e~xdx  = (xne~x/n) {1  + x/(n  + 1)  + x2/(n  + 1) (n  + 2) 

n 

+ x3/ (n  + 1)  (n  + 2)  (n  + 3)  + . 

This  series  also  converges  for  all  values  of  x,  but  is 
quicker  calculated  for  moderate  values  of  x than  the  pre- 
ceding series,  especially  if  n is  large. 

(c)  Again,  by  integrating  by  parts  show  that 


fxxn~1e’~x dx  = r(n)  - xn~1e~x{l  + (n  - l)/x  + (n-1) (n-2)/x2 
o 

+ (n-l) (n-2) (n-3)/x®  + ...  }. 

If  n is  an  integer,  this  series  will  terminate.  If  n is 
not  an  integer,  the  series  will  be  infinite  and  moreover 
will  diverge  for  any  x no  matter  how  large.  This  can  be 
seen  from  the  ratio  test,  since  the  ratio  of  the  s-th  term 
to  (s-l)th  term  is  (n-s)/x.  Since  n > 0,  the  series  will 
eventually  oscillate,  and  oscillate  infinitely  for  any  x. 
But  no  matter  where  curtailed,  the  error  committed  will  be 
less  than  the  next  term  . (Laplace  Theorie  analytique,  pp. 
174-5  in  the  3d  ed.,  1820).  There  will  always  be  an  opti- 
mum term  (Bayes,  1763),  and  the  series  carried  just  short 
of  the  optimum  term  will  give  the  best  possible  representa- 
tion for  those  values  of  x and  n.  When  x is  large,  a very 
few  terms  will  suffice  for  accurate  calculations,  especial- 
ly if  n is  not  large. 

Note : Series  (b)  and  (c)  can  also  be  derived 

by  the  D operator  used  inversely,  where 
D = d/dx  and  l/D  = /. 
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28.  Show  that  f xne~x  dx  = xne~xv,  where 
0 

v = 1 

1 - ny  1 

l + y i_ 

(y  = l/x)  1 - (n-l)y  * 

1 + 2y  £ 

1 - (n-2)y  1 

1 + 5y  i 

1+-. . . 

Hint: 

1 ' Show  that  v satisfies  the  differential  equation 

x dv/dx  = (x-n)v  - x.  (l) 

2 ’ Show  that  a solution  of  the  general  Riccati  equation 

x dv/dx  = (x-a)v  - x + bv2  (2) 

is  v = 1/(1  + kyvj  = x/(x  + kVj. ) , (y  = l/x)  (j>) 

provided  vx  satisfies  the  differential  equation 

x dVi/dx  = (x  + a + ljv!  - x (b-a)/k + kv^ . (4) 

3'  This  is  a Riocati  equation  like  (2)  if  k = b - a. 

Hence  a solution  of  (2)  is 

v = 1/(1  + [b-a] vxy)  (5) 

wherein  Vi  is  such  that  it  satisfies  the  differential 
equation  • 

x dVi/dx  = (x  + a+l)  - x + (b-a)vi.  (6) 

4'  It  is  evident,  then,  that 


Vi  = 1/(1  + [b  + l]yv2) 


(7) 
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if  v2  satisfies 

dv2/dx  = (x  - a)v2-x  + (b  + l)v22  (8) 

This  is  like  Eq . ( 2 ),  hut  with  h replaced  by  b + 1.  Hence 
we  may  say  that 


V = 

1/(1 

+ 

[b 

- 

a] 

yvi) 

Vi  = 

1/(1 

+ 

[b 

+ 

1] 

yv2) 

v2  = 

1/(1 

+ 

[b 

+ 

1 

- a]yv3) 

V3  = 

1/(1 

+ 

[b 

+ 

2] 

yv4) 

v4  = 

1/(1 

+ 

[b 

+ 

2 

- a]yv5) 

etc.,  and  the  required  result  is  established. 

29.  The  mean  of  the  curve  y = xn“^  e ~x  (Fig.  1, 

p.  2)  between  0 and  « is  at  n,  its  mode  (maximum)  is  at 

n - 1,  its  second  moment  coefficient  about  the  y axis  is 
n(n  + l),  and  its  standard  deviation  is  n?. 

30.  The  mean  of  the  curve  y = - x)n“ 1 lying 

between  0 and  1 is  at  m/(m  + n),  its  mode  is  atx(m  - l)/ 

(m  + n - 2),  and  its  standard  deviation  is  (mn)V 

(m  + n) (m  + n + 1)2. 


curve 


31.  Show  by  Eq.  (19)  that  if  the  area  under  the 

__  i 

ydz  = C(l  + z2)  2"  dz 


from  -oc  to  +oc  is  unity,  i.e.,  if  the  curve  is  "normalized," 
C must  have  the  value  l/B[J(n  - l),J]  . This  curve  is  the 
distribution  of  Student’s  z,  n being  the  sample  size. 


32.  Prove  that  the  standard  deviation  of  the  dis- 
tribution of  Student’s  z is  1/V(n  - J>) . (Student's  distri- 
bution of  z is  the  curve  used  in  the  preceding  exercise.) 
What  happens  to  the  standard  deviation  of  this  curve  if  n 
is  3 or  less? 
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33.  Helmert's  curve  for  the  distribution  of  the 
standard  deviation  s in  samples  of  n drawn  from  a normal 
universe  is 

y ds  = C (f)"-* 2  e'nsS/za2  ds 

If  the  area  under  the  curve  from  s = 0 to  s = oo  is  unity, 

ni(n-l) 

the  value  of  C must  be  — — , which  is  the 

r[i(n-l)]  2?(”-5)  a 

normalizing  factor.  The  distribution  of  s was  first  de- 
rived by  Helmert  in  1876;  see  p.  126  In  Deming  and  Birge's 
Statistical  Theory  of  Errors,  published  by  the  Graduate 

School,  Department  of  Agriculture. 

Hint:  In  this  exercise,  and  the  next  nine  it  will 

be  found  convenient  to  use  the  change  of  variable 
ns2/2 <y2  = x and  2 ds/s  = dx/x,  and  to  write  Eq. 

(2)  in  the  form 


34.  Show  that  if  P(s)  is  defined  by  the  equation 


P( 


L (|) 


n-2  _n  s2 


ds 


,»  ,S,n~2  -nsS/2C2 

fo  (?)  e ds 


then 


P(s)  = 1 - 


rjHn-1)] 
r ti(n-l)] 


The  integral  P(s)  is  the  area  under  the  distribution  of  s 
from  the  abscissa  s to  infinity,  and  the  evaluation  of 
p(s)  is  here  seen  to  depend  on  an  incomplete  Gamma  func- 
tion. See  Exercise  42  for  the  identification  of  P(s)  with 
p(x)  which  is  the  corresponding  integral  under  the  distri- 
bution of  X* 
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35.  Show  that  the  mean  of  the  curve  In  the  pre- 
ceding exercise  is 

1 

Bti(n-l),  J] 

The  symbol  E denotes  expectation  or  mathematical  average. 

Es  is  the  theoretical  average  standard  deviation  in  samples 
of  n drawn  from  a normal  universe  of  standard  deviation  a. 

36.  Show  that  for  Helmert's  distribution  of  s 
(Exercise  55) , 

Es2  = /°°  s2  y ds  = — a2 
0 n 


Es2  is  the  theoretical  average  of  the  square  of  the  stand- 
ard deviation  in  samples  of  n drawn  from  a normal  universe 
or  any  other  universe  of  standard  deviation  <7  . 


37.  Helmert's  curve  for  the  distribution  of  the 
root-mean-square  error  6 in  samples  of  n drawn  from  a 
normal  universe  is 


y d6  c 


-nh2/2GZ 

e 


d6 


If  the  area  under  the  curve  from  5 = 0 to  8 = ® is  unity. 


1 

,2« 


the  value  of  C must  be 
normalizing  factor. 


r(£n)  2=("-2)  a’ 


which  is  the 


38.  For  the  curve  in  the  preceding  exercise  show 


that 


E 8 " C 87  d8  - * bICTT 

I 52  = /"  82  J ® = cS 
0 

ES  is  the  theoretical  average  root-mean-square  error  In 
samples  of  n drawn  from  a normal  universe,  and  E52  is  the 
theoretical  average  mean  square  error  for  a normal  universe 
or  any  other  universe. 
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39.  Without  actually  performing  the  integrations 
or  reducing  the  integrals  to  Gamma  functions,  but  rather 
by  making  use  of  the  results  in  Exercises  36  and  38,  show 

that  the  mean  of  the  distribution  of  s2  must  be  — — O' 2 and 

n 

that  the  mean  of  the  distribution  of  8 2 must  be  a2. 


40.  When  there  are  k degrees  of  freedom  the  dis- 


tribution of  x 


y ax  = c (x2) 


SMk-2)  -if 


dx* 


If  the  area  under  the  curve  from  x = 0 to  x = « Is  to  be 


unity,  show  that  C must  have  the.  value 


r(Jk)  22* 


41.  The  mean  of  the  distribution  of  x Is  at  k and 
its  standard  deviation  is 

42.  With  P(s)  defined  as  in  Exercise  33  and  with 


P(X)  = T7  S 

r(£k)  22k  x 


00  / 2 \ 2 ( k-2 ) ■ - 1 X 


1 „,2 


(x2) 


dx‘ 


show  that  if  ns2/a2  = x2  and  n - 1 = k,  then  P(s)  = P(x). 
The  results  of  this  example  show  that  instead  of  finding 
P(s)  from  the  Tables  of  the  Incomplete  Gamma  Function, 

P(s)  can  be  found  from  tables  of  which  give  P(x)  as  a 
function  of  X2  and  k. 

43.  Show  that  the  value  of  C required  to  normalize 
the  curve 


, n f3\n-l  -ns2/ 2<?2  * 
y da  = C (-)  e d<j 


i(n-2) 


from  a = 0 to  a = co  is 


r[i(n-2) ] 2=("-4)  s 


44.  Show  that  the  value  of  C required  to  normalize 
the  curve 

, „ ( a x c+  5 - o2/  202  , _ 

y dc  = C (jj)  © d(T 
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1 

from  a=Otoa=ooig1  . (This  curve  is 

2*c  a r (ic  + 1) 

Molina  and  Wilkinson* s suggestion  for  a "prior  existence" 
curve  of  ff ; Bell  System  Technical  Journal,  vol.  8,  1929: 
pp.  652-645.) 

45 . The  mean  of  the  curve  in  the  preceding  exercise 
is  at  Eff  = [si/y(2n)]  B(J  c+1,  J),  and  its  second  moment 
coefficient  about  the  origin  is  Ecr2  = a2/c. 

46 . The  curve 

1-tt2 

y = yo(l+at)  a e*t/<a 
a ~2 

Yo  = a/(a2e)  r(a"2) 

is  sometimes  used  to  approximate  the  point  binomial.  This 
is  equivalent  to  a Pearson  Type  III  curve.  A table  of  in- 
tegrals for  this  curve  was  published  by  Salvosa.13 

(a)  Verify  that  the  value  given  for  y0  actually  makes  it  a 
normalizing  factor. 

(b)  Show  that  the  mean  of  this  curve  is  at  t = 0;  in  other 
words,  the  origin  is  taken  at  the  mean. 

(c)  The  second  moment,  coefficient  of  the  curve  about  the 
origin  is  unity;  hence  the  S.D.  is  unity. 

(d)  The  third  moment  coefficient  of  the  curve  is  2a. 

(e)  The  mode  (maximum)  is  at  -a. 

/ \ mean— mode 

(f)  If  skewness  be  defined  as  — - — — — -,  then  the  skewness 

S . D • 

of  this  curve  is  a. 

(g)  If  skewness  be  defined  as  n3/2(S.D. )3,  then  the  skew- 
ness of  this  curve  is  a.  ji 3 denotes  the  third  moment 
coefficient  of  the  curve  about  its  mean.  The  two 
measures  of  skewness  just  given  happen  to  coincide  for 
the  Type  III  curve,  but  for  no  other  curve  with  skew- 
ness different  from  zero. 


